Numerical study of hot strongly interacting matter 

1 P. Petreczky 

Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 

Abstract. I review recent progress in study of strongly interacting matter at high 
temperatures using Monte-Carlo simulations in lattice QCD. 

^"'l 1. Introduction 

It is expected that strongly interacting matter undergoes a transition in some temperature 
interval from hadron gas to deconfined state also called the quark gluon plasma (QGP) [lj. 
Creating deconfined medium in a laboratory is the subject of the large experimental program at 
RHIC [2] and LHC [3]. Attempts to study QCD thermodynamics on the lattice go back to the 
early 80's when lattice calculations in SU(2) gauge theory provided the first rigorous theoretical 
t^j- ' evidence for deconfinement [U El E] • The problem of calculating thermodynamic observables in 

pure gluonic theory was solved in 1996 [7], while simulations involving dynamical quarks were 
t^J" ! limited to large quark masses and had no control over the discretization errors [SJ El ED] . Lattice 

calculations of QCD thermodynamics with light dynamical quarks remained challenging until 
recently. During the past six years calculations with light u, d quarks have been performed using 
improved staggered fermion actions [HI HH H31 [TH [151 Ull UHl [HI EQl EU E21 El]. There 
was also progress in lattice QCD calculations at non-zero temperature using Domain Wall and 
^ , Wilson fermion formulations |25 } 126 } 127], However, due to much larger computational costs the 

corresponding results are far less extensive. 

To get reliable predictions from lattice QCD the lattice spacing a should be sufficiently small 
relative to the typical QCD scale, i.e. Kqcdcl <C 1. For staggered fermions, discretization errors 
go like 0((ciAqcd) 2 ) but discretization errors due to flavor symmetry breaking turn out to be 
numerically quite large and dominate the cutoff dependence of thermodynamic quantities at low 
temperatures. To reduce these errors one has to use improved staggered fermion actions with 
so-called fat links [28]. At high temperature the dominant discretization errors go like (aT) 2 
and therefore could be very large. Thus it is mandatory to use improved discretization schemes, 
which improve the quark dispersion relation and eliminate these discretization errors. Lattice 
fermion actions used in numerical calculations typically implement some version of fat links as 
well as improvement of quark dispersion relation and are referred to as p4, asqtad, HISQ /tree 
and stout. 

In this contribution I am going to discuss lattice QCD calculations on the equation of state, 
study of deconfinement aspects of the QCD transition, including color screening and fluctuations 
of conserved charges, determination of the chiral transition temperature, as well as calculations 
of temporal and spatial meson correlation functions. 
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2. Equation of State 

The equation of state has been calculated using different improved staggered fermion actions 
pA, asqtad, stout and HISQ/tree. In the lattice calculations of the equation of state and 
many other quantities the temperature is varied by varying the lattice spacing at fixed value of 
the temporal extent N T . The temperature T is given by the lattice spacing and the temporal 
extent, T = l/{N T a). Therefore taking the continuum limit corresponds to iV T — > oo at the 
fixed physical volume. The calculation of thermodynamic observables proceeds through the 
calculation of the trace of the energy momentum tensor, e — 3p, also known as trace anomaly 
or interaction measure. This is due to the fact that this quantity can be expressed in terms 
of expectation values of local gluonic and fermionic operators, (see e.g. Ref. [IS]). Different 
thermodynamic observables can be obtained from the interaction measure through integration 
of the trace anomaly 0. The pressure can be written as 

p(t) p(t ) fT d r 

The lower integration limit To is chosen such that the pressure is exponentially small there. 
Furthermore, the entropy density can be written as s = (e +p)/T. Since the interaction 
measure is the basic thermodynamic observable in the lattice calculations it is worth to discuss 
its properties more in detail. In Fig. Q] (left panel) I show the results of calculation with pA 
and asqtad actions using N T = 6 and 8 lattices and light quark masses mi = m s /10, where m s 
is the physical strange quark mass. These calculations correspond in the continuum limit to 
the pion mass of 220MeV and 260MeV for pA and asqtad respectively. The interaction measure 
shows a rapid rise in the transition region and after reaching a peak at temperatures of about 
200MeV decreases. Cutoff effects (i.e. N T dependence) appears to be the strongest around the 
peak region and decrease at high temperatures. For temperatures T < 270MeV calculations 
with pA and asqtad actions have been extended to smaller quark masses, mi = m s /20, that 
correspond to the pion mass of about 160MeV in the continuum limit |171 129] . It turns out 
that the quark mass dependence is negligible for mi < m s /10. Furthermore, for astqad action 
calculations have been extended to N T = 12 lattices [29]. The trace anomaly was calculated 
with HISQ/tree action on lattices with temporal extent N T = 6 and 8 and mi = m s /20 [29] 
(corresponding to = 160MeV in the continuum limit). Finally, calculation of the trace 
anomaly and the equation of state was performed with stout action using N T = 4, 6, 8, 10 and 
12 and physical light quark masses [21]. Using the lattice data from N T = 6, 8, 10 a continuum 
estimate for different quantities was given [21] . In Fig. [Tj (right panel) the results of different 
lattice calculations of e — 3p are summarized corresponding to the pion masses close to the 
physical value. I also compare the lattice results with the parametrization s95p-vl of e — 3p. This 
parametrization combines lattice QCD results of Refs. |16tll8j at high temperatures with hadron 
resonance gas model (HRG) at low temperatures (T < 170MeV) [31] . At low temperatures there 
is a fair agreement between the results obtained with stout action and the results obtained with 
HISQ/tree action as well as with asqtad action for N T = 12. All these lattice results are slightly 
above the HRG curve. There is a big difference in the height of the peak in the calculations 
obtained with asqtad and HISQ/tree actions and the continuum estimate for stout action. The 
pA action has large cutoff effects at low temperatures and in the peak region. As the result the 
pA results are higher in the peak region and fall below the HRG curve at low temperatures. 
When cutoff effects in the hadron spectrum are taken into account in the HRG model a good 
agreement between the pA data and HRG result can be achieved [31]. Since the dominant cutoff 
effects of order a 2 T 2 are eliminated, the iVV-dependence is expected to be small for pA, asqtad 
and HISQ/tree action at high temperatures. We see that for T > 250MeV all these lattice 



1 A somewhat different approach was used in Ref. [21] 



actions lead to similar results. At temperatures above 350MeV we also see a good agreement 
with the stout results. 
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Figure 1. The interaction measure calculated with mi = m s /10 and pA and asqtad actions 
|18j (left) and with mi = m s /20 for pA action [17] as well as with HISQ/tree and asqtad actions 
|29j . Also shown in the figure are the continuum estimates obtained with stout action and the 
parametrization based on hadron resonance gas (HRG) model |31j . 

The pressure, the energy density and the entropy density are shown in Fig. [2j The energy 
density shows a rapid rise in the temperature region (170 — 200) MeV and quickly approaches 
about 90% of the ideal gas value. The pressure rises less rapidly but at the highest temperature 
it is also only about 15% below the ideal gas value. In the previous calculations with the pA 
action it was found that the pressure and the energy density are below the ideal gas value by 
about 25% at high temperatures [10j . Possible reason for this larger deviation could be the fact 
that the quark masses used in this calculation were fixed in units of temperature instead being 
tuned to give constant meson masses as lattice spacing is decreased. As discussed in Ref. [30] 
this could reduce the pressure by 10 — 15% at high temperatures. In Fig. [2] I also show the 
entropy density divided by the corresponding ideal gas value and compare the results of lattice 
calculations with resummed perturbative calculation [32, 33] as well as with the predictions from 
AdS/CFT correspondence for the strongly coupled regime |34j . The later is considerably below 
the lattice results. Note that the pressure, the energy density and the trace anomaly have also 
been recently discussed in the framework of resummed perturbative calculations which seem to 
agree with lattice data quite well at high temperatures [35 . 

The differences between the stout action and the pA and asqtad actions for the trace anomaly 
translates into the differences in the pressure and the energy density. In particular, the energy 
density is about 20% below the ideal gas limit for the stout action. 

3. Deconfinement : Fluctuations of conserved charges 

Due to the infamous sign problem lattice QCD Monte-Carlo simulations are not possible at non- 
zero quark chemical potentials. The pressure at non-zero chemical potentials can be evaluated 
using Taylor expansion. The Taylor expansion can be set up in terms of quark chemical potentials 
or in terms of chemical potentials corresponding to baryon number B, electric charge Q and 
strangeness S of hadrons. The expansion coefficients in quark chemical potential x^uls anc ^ 
the hadronic ones x'bqs are re l a ted to each other [36]. While Taylor expansion can be used 
to study the physics at non-zero baryon density they are interesting on their own right as 
they are related to the fluctuations of conserved charges. The later are sensitive probes of 
deconfinement. This is because fluctuation of conserved charges are sensitive to the underlying 
degrees of freedom which could be hadronic or partonic. Fluctuations of conserved charges have 



0.4 0.6 0.8 1 1.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 




AdS/CFT 



100 150 200 250 300 350 400 450 500 550 100 200 300 400 500 600 700 800 

Figure 2. The energy density and the pressure as function of the temperature (left), and the 
entropy density divided by the corresponding ideal gas value (right). The dashed lines in the 
right panel correspond to the resummed perturbative calculations while the solid black line is 
the AdS/CFT result. 



been studied using different staggered actions [IT] [18] [20], [36], [37]. As an example in Fig. [3] I 
show the quadratic strangeness fluctuation as the function of the temperature. Fluctuations 
are suppressed at low temperatures because conserved charges are carried by massive strange 
hadrons (mostly kaons). They are well described by Hadron Resonance Gas (HRG) model at low 
temperatures. Strangeness fluctuations rapidly grow in the transition region of T = (160 — 200) 
MeV as consequence of deconfinement. At temperatures T > 250 MeV strangeness fluctuations 
are close to unity that corresponds to the massless ideal quark gas. Similar picture can seen in 
other fluctuations, in particular for light quark number fluctuation, Xh which is also show in 
Fig. The more rapid rise of Xl m the transition region is presumably a quark mass effect. 

The transition from hadronic to quark degrees of freedom can be particularly well seen in the 
temperature dependence of the Kurtosis, which is the ratio of quartic to quadratic fluctuations. 
This quantity can be also measured experimentally. In Fig. [5] I show the Kurtosis of the baryon 
number as a function of the temperature. At low temperatures it is temperature independent 
and is close to unity in agreement with the prediction of the HRG model. In the transition 
region it sharply drops from the hadron resonance gas value to the value corresponding to an 
ideal gas of quarks. Since the fluctuations of conserved charges are so well described by ideal 
quark gas it is interesting to compare the lattice results for quark number fluctuations with 
resummed perturbation theory. Quadratic quark number fluctuations are also called quark 
number susceptibilities and have been calculated in resummed perturbation theory [38 ] 139 ] 14*0]. 
The comparison of the resummed perturbative results with lattice calculations performed with 
pA action [UJ for the quark number susceptibility divided by the ideal gas value is shown in 
Fig. HI The figure shows a very good agreement between the lattice calculations and resummed 
perturbative results. Let me finally note that the discretization errors in the lattice calculations 
of the fluctuations are small at high temperatures when calculated with the p4 action. The same 
is true for asqtad action [TT] [18] . 

Fluctuations of conserved charges have also been calculated using gauge/string duality in 
Refs. [42] 143] . These calculations could reproduce several features of the lattice data. 

4. Deconfinement : color screening 

One of the most prominent feature of the quark gluon plasma is the presence of chromoelectric 
(Debye) screening. The easiest way to study chromoelectric screening is to calculate the Polyakov 
loop. The Polyakov loop is an order parameter for the deconfinement transition in pure gauge 
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Figure 3. Strangeness fluctuation calculated with asqtad and HISQ/tree actions |24j and 
compared with the continuum estimate obtained with stout action |20j as well as the prediction 
of the HRG model shown as a black line (left panel). In the right panel the light quark number 
fluctuations are shown for HISQ/tree action |24| . 
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Figure 4. The Kurtosis of baryon number (left) and the quark susceptibility compared with 
resummed perturbative calculations (right). The open symbols on the right plot correspond to 
asqtad results for strange quark number susceptibility [TT]. The resummed perturbative results 
were obtained in next-to- leading approximation (NLA) |39j . 

theory, which is governed by Z(N) symmetry. For QCD this symmetry is explicitly broken 
by dynamical quarks. There is no obvious reason for the Polyakov loop to be sensitive to the 
singular behavior close to the chiral limit although speculations along these lines have been 
made [52]. The Polyakov loop is related to the screening properties of the medium and thus 
to deconfinement. After proper renormalization, the square of the Polyakov loop characterizes 
the long distance behavior of the static quark anti-quark free energy; it gives the excess in free 
energy needed to screen two well-separated color charges. The renormalized Polyakov loop has 
been studied in the past in pure gauge theory [5IJ [53] as well as in QCD with two [51], three 
[50] and two plus one flavors [Ml HE]. The renormalized Polyakov loop, calculated on lattices 
with temporal extent N T , is obtained from the bare Polyakov 



L ren (T) = z((3) N ^ L bare (P) = z(f3) 



^TrW(x) 



N T -1 

W(x)= J] U (xo,x), 

x =0 



(2) 



where Uq = exp(igaAo) denotes the temporal gauge link and z{j3) is the renormalization constant 
determined from the T = static potential [24J. The numerical results for the renormalized 
Polyakov loop for the HISQ/tree action are shown in Fig. [5J As one can see from the figure the 
cutoff (N T ) dependence of the renormalized Polyakov loop is small. We also compare our results 
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Figure 5. The Polyakov loop as function of the temperature in pure gauge theory and 2+1 
flavor QCD. 



with the continuum extrapolated stout results [20] and the corresponding results in pure gauge 
theory [511 [53]. We find good agreement between our results and the stout results. We also 
see that in the vicinity of the transition temperature the behavior of the renormalized Polyakov 
loop in QCD and in the pure gauge theory is quite different. 

Further insight on chromoelectric screening can be gained by studying the singlet free energy 
of static quark anti-quark pair (for reviews on this see Ref. |55[ [56] ) , which is expressed in terms 
of correlation function of temporal Wilson lines in Coulomb gauge 

exp(-F 1 (r,r)/T) = ^Tr(W(r)W^0)). (3) 

Instead of using the Coulomb gauge the singlet free energy can be defined in gauge invariant 
manner by inserting a spatial gauge connection between the two Wilson lines. Using such 
definition the singlet free energy has been calculated in SU(2) gauge theory [57J. It has been 
found that the singlet free energy calculated this way is close to the result obtained in Coulomb 
gauge [57]. The singlet free energy turned out to be useful to study quarkonia binding at high 
temperatures in potential models (see e.g. Refs. [HI H5[ SSJ, H7J 48J). The singlet free energy 
also appears naturally in the perturbative calculations of the Polyakov loop correlators at short 
distances [58] . 

The singlet free energy was recently calculated in QCD with one strange quark and two 
light quarks with masses corresponding to pion mass of 220MeV on 16 3 x 4 lattices 09]. The 
numerical results are shown in Fig. [6J At short distances the singlet free energy is temperature 
independent and coincides with the zero temperature potential. In purely gluonic theory the 
free energy grows linearly with the separation between the heavy quark and anti-quark in the 
confined phase. In presence of dynamical quarks the free energy is saturated at some finite value 
at distances of about 1 fm due to string breaking [551 1501 154"] . This is also seen in Fig. [6] Above 
the deconfinement temperature the singlet free energy is exponentially screened at sufficiently 
large distances |51[ 153] with the screening mass proportional to the temperature , i.e. 

Fi r,T = FooCT) - ^^exp(-m D (T)r), mD ~ T. (4) 
3 47rr 

Therefore in Fig. 01 also show the combination F\(r,T) — F 00 (T) as a function of rT. As one 
can see from the figure this function shows an exponential fall-off at distances rT > 0.8. The 
fact that the slope is the same for all temperatures means that mo ~ T, as expected. 
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Figure 6. The singlet free energy F\(r, T) calculated in Coulomb gauge on 16 3 x 4 lattices 
(left) and the combination F\(r,T) — Foq(T) as function of rT (right). The solid black line is 
the parametrization of the zero temperature potential. 



Let me finally note that contrary to the electro magnetic plasma the static chormomagnetic 
fields are screened in QGP. This is due to the fact that unlike photons gluons interact with each 
other (the stress tensor is non-linear in QCD). Magnetic screening is non-perturbative, i.e. it 
does not appear at any finite order of pertubation theory. In lattice calculations chromomagnetic 
screening is studied either in terms of spatial Wilson loops [59] or in terms of spatial gluon 
propagators |60t [6T1 I62j . The numerical results obtained so far show that the length scale 
related to magnetic screening is larger than the one related to electric screening. 

5. Chrial transition 

The Lagrangian of QCD has an approximate SUa(3) chiral symmetry. This symmetry is broken 
in the vacuum. The chiral symmetry breaking is signaled by non-zero expectation value of 
the quark or chiral condensate, (ijjip) ^ in the massless limit. This symmetry is expected 
to be restored at high temperatures and the quark condensate vanishes. There is an explicit 
breaking of the chiral symmetry by the small value of u, d and s quark masses. While due to the 
relatively large strange quark mass (m s ~ 100 MeV) SUa(3) may not be a very good symmetry 
its subgroup SUa(2) remains a very good symmetry and is relevant for the discussion of the 
finite temperature transition in QCD. If the relevant symmetry is SUa(2) the chiral transition is 
expected to be second order for massless light (u and d) quarks belonging to the 0(4) universality 
class. Recent calculations with p4 action support this pictures [63]. This also means that for non- 
zero light quark masses the transition must be a crossover. The later fact seems to be supported 
by calculations in Ref. [64] • The Ua(1) symmetry is explicitly broken in the vacuum by the 
anomaly but it is expected to be effectively restored at high temperatures as non-perturbative 
vacuum fluctuations responsible for its breaking are suppressed at high temperatures. If the 
Ua(X) symmetry is restored at the same temperature as the SUa(2) symmetry the transition 
could be first order [65] . Recent calculations with staggered as well as with domain wall fermions 
suggest that J7a(1) symmetry gets restored at temperature that is significantly higher than the 
chiral transition temperature [25\ 166]. 

For massless quark the chiral condensate vanishes at the critical temperature T® and is the 
order parameter. Therefore in the lattice studies one calculates the chiral condensate and its 
derivative with respect to the quark mass called the chiral susceptibility. For the staggered 
fermion formulation most commonly used in the lattice calculations these quantities can be 



written as follows: 



(rp\ d(lf;ip) qiT 

Xm,q{T) = Uf = Xq,disc + Xq,con ? = M> (6) 

where the subscript x = r and x = will denote the expectation value at finite and zero 
temperature, respectively. Furthermore, D q = m q ■ 1 + D is the fermion matrix in the canonical 
normalization and nj = 2 and 1 for light and strange quark. In Eq. ([6] we made explicit 
that chiral susceptibility is the sum of connected and disconnected Feynman diagrams. The 
disconnected and connected contributions can be written as 



Xq,di 



"7 



{((TrD- 1 ) 2 ) - (TrD- 1 ) 2 } , (7) 
Xq,con = -^Tt^D-^x^D-^x)) , q = l,s. (8) 



16N*N 7 
4 



The disconnected part of the light quark susceptibility describes the fluctuations in the light 
quark condensate and is directly analogous to the fluctuations in the order parameter of an 
O(N) spin model. The second term (xq,con) arises from the explicit quark mass dependence 
of the chiral condensate and is the expectation value of the volume integral of the correlation 
function of the (isovector) scalar operator tpip. Let me note that in the massless linit only %i disc 
diverges. 



5.1. The temperature dependence of the chiral condensate 

The chiral condensate needs a multiplicative, and also an additive renormalization if the quark 
mass is non-zero. Therefore the subtracted chiral condensate is considered 

lA j ~<^,o-^W> s ,o' (Jj 

In Fig. [7] I show the results for Aj s calculated with HISQ/tree action and compared to the 
renormalized Polyakov loop and light quark number fluctuation discussed in relation to the 
deconfining transition. Interestingly the rapid decrease in the subtracted chiral condensate 
happens at temperatures that are smaller than the temperatures where the Polyakov loop 
rises rapidly. On the other hand the rapid change in \i an d A/ s happen roughly in the same 
temperature region. 

Another way to get rid of the multiplicative and additive renormalization is to subtract the 
zero temperature condensate and multiply the difference by the strange quark mass, i.e. consider 
the following quantity 

Af = d + n f m s rf((^) q , T - (^) q , ) q = l,s. (10) 

As before, nf = 2 for light quarks and rif = 1 for strange quarks, while d is a normalization 
constant. The factor rf was introduce to make the combination dimensionless. Here r\ is the 
scale parameter defined from the zero temperature static potential [24J. It is convenient to 
choose the normalization constant to be the light quark condensate for mi = multiplied by 
m s rf. In Fig. [8] the renormalized quark condensate is shown as function of temperature for 
HISQ/tree and stout actions. We see a crossover behavior for temperature of (150 — 160)MeV, 
where A^ drops by 50%. The difference between the stout and HISQ/tree results is a quark 
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Figure 7. The subtracted chiral condensate calculated with HISQ/tree action on N T = 8 
lattice and compared to the renormalized Polyakov loop (left) and light quark number fluctuation 
(right). 



0.025 



0.015 
0.01 



- a 


HISQ/tree : N x =1 2 ■ 




N T =8 □ 


A l * 


N T =6 x 


N T =8, m|=0.037m s ♦ 
stout, cont. •» 


rax 


*»> 




#♦ m 












**i 


" i 


T 


[MeV] 



120 140 160 180 200 220 240 



0.025 
0.02 

0.015 
0.01 

0.005 


-0.005 



Si 



m HISQ/tree: N T =1 2 ■ 
% N T =8 

B* N.=6 x 

■ E 



T [MeV] 

120 140 160 180 200 220 240 



Figure 8. The renormalized chiral condensate for the HISQ/tree action with mi/m s = 0.05 
is compared to the stout data. In the right panel, we show the renormalized strange quark 
condensate for the HISQ/tree action. 



mass effect. Calculations for HISQ/tree action were performed for m n = 160 MeV, while the 
stout calculations were done for the physical quark mass. For a direct comparison with stout 
results, we extrapolate the HISQ/tree data in the light quark mass and also take care of the 
residual cutoff dependence in the HISQ/tree data. This was done in Ref. |24] and the results 
are shown in the figure as black diamonds demonstrating a good agreement between HISQ/tree 
and stout results. Contrary to A^ the renormalized strange quark condensate A^ shows only a 
gradual decrease over a wide temperature interval dropping by 50% only at significantly higher 
temperatures of about 190 MeV. 



5.2. The crial susceptibility 

For a true chiral phase transition the chiral susceptibility diverges at the transition temperature. 
For physical value of the quark masses we expect to see a peak in the chiral susceptibility at 
certain temperature that defines the crossover temperature. The chiral susceptibility also needs 
a multiplicative and additive renormalization Therefore the following quantity is considered 

^jP- = (Xm,l(T) - Xm,l(T = 0)) . (11) 



The numerical results for this quantity were presented in Ref. |14j for stout action and in Ref. 
|24j for the HISQ /tree action. It was shown that the results of the two calculations agree quite 
well up to small quark mass effects [24J 0. 

5.3. 0(N) scaling and the transition temperature 

In the vicinity of the chiral phase transition, the free energy density may be expressed as a sum 
of a singular and a regular part, 



T 

f = ~y In Z = f sing {t, h) + f reg {T, m u m s ) . (12) 



Here t and h are dimensionless couplings that control deviations from criticality. They are 
related to the temperature T and the light quark mass mi as 

<=i^ . H = ±.H , H=?L, (13) 
t T c u h m s 

where T® denotes the chiral phase transition temperature, i.e., the transition temperature at 
H = 0. The scaling variables t, h are normalized by two parameters to and ho, which are unique 
to QCD and similar to the low energy constants in the chiral Lagrangian. These need to be 
determined together with T°. In the continuum limit, all three parameters are uniquely defined, 
but depend on the value of the strange quark mass. 

The singular contribution to the free energy density is a homogeneous function of the two 
variables t and h. Its invariance under scale transformations can be used to express it in terms 
of a single scaling variable 

i/ps _ 1 T - T c (ho\ 1/P5 1 T ~ T r ( 1 \ 17/35 



where /3 and 5 are the critical exponents of the O(N) universality class and zo = to/h^j^ 5 . Thus, 
the dimensionless free energy density / = f /T 4 can be written as 

f(T,m h m s ) = h 1+1 / 5 f f (z) + f reg (T,H,m s ) , (15) 

where ff is the universal scaling function and the regular term f reg gives rise to scaling violations. 
This regular term can be expanded in a Taylor series around (£, h) = (0, 0). 

It should be noted that the reduced temperature t may depend on other couplings in the 
QCD Lagrangian which do not explicitly break chiral symmetry. In particular, it depends on 
light and strange quark chemical potentials (J, q , which in leading order enter only quadratically, 



+ Ki s — — . (16) 





The transition temperature can be defined as peaks in susceptibilities (response functions) 
that are second derivatives of the free energy density with respect to relevant parameters. Since 
there are two relevant parameters we can define three susceptibilities: 

d 2 f d 2 f 8 2 f 

Xm > l -dmf> Xt ' l ~dW^i' Xt *-W (17) 



2 In Ref. [14] the light quark mass was used instead of m B in Eq. (jlll) . In the comparison with was taken into 
account [24| . 



Thus three different pseudo-critical temperatures T m j, T t j and i can be defined. In the vicinity 
of the critical point the behavior of these susceptibilities is controlled by three universal scaling 
function that can be derived from ff. In the chiral limit T m j = T t j = T t t = T®. There is, 
however, an additional complication for O(N) universality class: while Xm,i and Xt,l diverge at 
the critical point for mi — > 

Xm,i^rn x ' , xt,l~m\ , (18) 

Xt.t is finite because a < for O(N) models ( xt,t ~ )• Therefore, one has to consider the 
third derivative of / with respect to t : 

d 3 f 

Xt,t,t = (19) 

In the vicinity of the critical point the derivatives with respect to t can be estimated by 
taking the derivatives with respect to {if, i.e. the response functions Xt,l and Xt,t,t are identical 
to the first Taylor expansion coefficient of the quark condensate and the sixth order expansion 
coefficient to the pressure, respectively. The former controls the curvature of the transition 
temperature as function of the quark chemical potential n q and was studied for p4 action using 
N T = 4 and 8 lattices [67] • The later corresponds to the sixth order quark number fluctuation 
which is related to the deconfinement aspects of the transition. The fact that this quantity is 
sensitive to the chiral dynamics points to a relation between deconfining and chiral aspects of 
the transition. In the following I discuss the determination of the transition temperature defined 
as peak position of Xm,U i-e. T c = T m j. 

5.4- Determination of the transition temperature 

The O(N) scaling described in the above subsection can be used to determine the pseudo- 
critical temperature of the chiral transition. For the study of the 0(N) scaling it is convenient 
to consider the dimensionless order parameter 

M b = m s ^^. (20) 

The subscript " b" refers to the fact that this is a bare quantity since the additive UV divergence 
is not removed. From the point of view of the scaling analysis this divergent term is just a regular 
contribution. For sufficiently small quark mass and in the vicinity of the transition region we 
can write 

M b (T, H) = hV'fait/hW 5 ) + f M ,reg(T, H). (21) 

Here fc{ z ) is the scaling function related to ff and was calculated for 0(2) and 0(4) spin models 
|68} [69l [70] . The regular contribution can be parametrized as 

f M ,reg(T,H) = a t {T)H 

a + a 1 ^j^- + a 2 [^^) | H. (22) 

Then we have the following behavior for the light chiral susceptibility 

Xm,l _ T 2 ( 1 ,1/5-1 f ( \ , df M ,reg{T,H) 

T2 " m^Uo h[) OH 

with f x (z) = ±[f G (z)-lf' G {z)]. (23) 




One then performs a simultaneous fit to the lattice data for Mf, and Xm,i treating T£, to, ho, ao, ai 
and a2 as fit parameters [23]. This gives a good description of the quark mass and temperature 
dependence of Xm,i and allows to determine accurately the peak position in Xm,l- As an example 
in Fig. I show the 0(4) scaling fits for N T = 8 lattice data obtained with HISQ/tree action. 
The scaling fit works quite well. Similar results have been obtained for N T = 6 and 12 as well as 
for asqtad action on N T = 8 and 12 lattices [24]. Furthermore, scaling fits have been performed 
assuming 0{2) universality class. The quality of these fits were similar to the 0(4) ones and 
the resulting transition temperatures turned out to be the same within statistical errors [24] . 
Having determined T c for HISQ/tree and asqtad action for each N T a combined continuum 
extrapolation was performed using different assumption about the N T dependence, of T c This 
analysis gave [23] : 

T c = (159 ± 9)MeV. (24) 

The analysis also demonstrated that HISQ/tree and asqtad action give consistent results 
in the continuum limit. The Budapest- Wuppertal collaboration found T c = 147(2)(3)MeV, 
157(3)(3)MeV and 155(3)(3)MeV defined as peak position in XRi inflection points in and 
respectively [20]. These agree with the above value within errors. 



6. Temporal meson correlars and spectral functions 

Information on hadron properties at finite temperature as well as the transport coefficients are 
encoded in different spectral functions. In particular the fate of different quarkonium states in 
QGP can studied by calculating the corresponding quarkonium spectral functions (for a recent 
review see Ref. [56] )• On the lattice we can calculate correlation function in Euclidean time. 
This is related to the spectral function via integral relation 

G(r,T) = r dw{u,,T)K{T, U ,T), K(r,u,T) = C ° sh( ^ ~- \^ )] . (25) 
Jo smh(uj/21 ) 

Given the data on the Euclidean meson correlator G(t, T) the meson spectral function can 
be calculated using the Maximum Entropy Method (MEM) [72]. For charmonium this was 
done by using correlators calculated on isotropic lattices [731 El] as wen as anisotropic lattices 
|75[ [76l [77] in the quenched approximation. It has been found that quarkonium correlation 




Figure 9. Scaling fits and data for the chiral condensate Mf, calculated with the HISQ/tree 
action on lattices with temporal extent N T = 8 (left) and the chiral susceptibility Xm,i (right). 
The data for Mf, at mi/m s = 0.025 and for Mb and Xm,i at mi/m s = 0.05 are fit simultaneously 
using the 0(4) scaling Ansatz. The points used in the scaling fits are plotted using open symbols. 
The dotted lines give the data scaled to the physical quark masses. 



function in Euclidean time show only very small temperature dependence [74 [ 177] . In other 
channels, namely the vector, scalar and axial-vector channels stronger temperature dependence 
was found [THE?]. The spectral functions in the pseudo-scalar and vector channels reconstructed 
from MEM show peak structures which may be interpreted as a ground state peak |74 [ 1751 176] . 
Together with the weak temperature dependence of the correlation functions this was taken 
as strong indication that the IS charmonia (r/ c and J/VO survive in the deconfined phase to 
temperatures as high as 1.6T C [7H [75j [76]. A detailed study of the systematic effects show, 
however, that the reconstruction of the charmonium spectral function is not reliable at high 
temperatures [77], in particular the presence of peaks corresponding to bound states cannot 
be reliably established. Presence of large cutoff effects at high frequencies also complicates the 
analysis [80J. The only statement that can be made is that the spectral function does not 
show significant changes within the errors of the calculations. Recently quarkonium spectral 
functions have been studied using potential models and lattice data for the singlet free energy 
of static quark anti-quark pair |46[ [4"7] 148] . These calculations show that all charmonium 
states are dissolved at temperatures smaller than 1.2T C , but the Euclidean correlators do not 
show significant changes and are in fairly good agreement with available lattice data both for 
charmonium |74[ [77j and bottomonium |77[ [75] . This is due to the fact that even in absence 
of bound states quarkonium spectral functions show significant enhancement in the threshold 
region [35]. Therefore previous statements about quarkonia survival at high temperatures have 
to be revisited. Exploratory calculations of the charmonium correlators and spectral functions 
in 2-flavor QCD have been reported in Ref. |79j and the qualitative behavior of the correlation 
functions was found to be similar. 

The large enhancement of the quarkonium correlators above deconfinement in the scalar and 
axial-vector channel can be understood in terms of the zero mode contribution |45[ IBT] and 
not due to the dissolution of the IP states as previously thought. Similar, though smaller in 
magnitude, enhancement of quarkonium correlators due to zero mode is seen also in the vector 
channel [77] . Here it is related to heavy quark transport |44[ [82] . Due to the heavy quark mass 
the Euclidean correlators for heavy quarkonium can be decomposed into a high and low energy 
part G(t, T) = G\ ow (t, T) + Ghigh(T, T) The area under the peak in the spectral functions at zero 
energy u ~ giving the zero mode contribution to the Euclidean correlator is proportional to 
some susceptibility, G 1 1ow (t,T) ~ Tx l (T), which have been calculated on the lattice in Ref. [83] . 
It is natural to ask whether the generalized susceptibilities can be described by a quasi-particle 
model. The generalized susceptibilities have been calculated in Ref. [84] in the free theory. 
Replacing the bare quark mass entering in the expression of the generalized susceptibilities by 
an effective temperature dependent masses one can describe the zero mode contribution very 
well in all channels [83]. 

While temporal correlators are not sensitive to the change in the spectral functions 
spatial quarkonium correlation functions could be more sensitive to this. Recent lattice 
calculations show indication for significant change in spatial charmonium correlators above 
deconfinement [85] . 

The analysis of the meson spectral functions was also performed in the light quark sector 
[86] [87[ l89l 155] . Somewhat surprisingly this analysis revealed peak structures in the deconfined 
phase which are difficult to interpret. Most likely the peak structures are artifacts of the MEM 
analysis. The new analysis of the vector spectral function performed on fine lattices and using 
several lattice spacings did not show any evidence for such structures [90]. Furthermore, the 
same analysis revealed the expected transport peak at low frequencies and estimated the electric 
conductivity to be [90] : 

1/3 < — ^- < 1, C em = e 2 Y,Ql (26) 



7. Conclusions 

Recently significant progress has been made in lattice QCD calculations at non-zero temperature. 
Chiral and deconfining aspects of the QCD transition have been studied using improved 
staggered quark formulation allowing to control discretization effects. It has been shown that 
in the continuum limit different discretization schemes give consistent results. In particular, 
agreement has been reached on the value of the chiral transition temperature. There is still 
disagreement in the lattice calculation of the equation of state. Finally significant progress has 
been made in lattice calculations of meson spectral functions that encode in-medium meson 
properties and transport coefficients. 
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